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The  Air  Force  General  Perturbation  theory  (AFGP4)  based 
on  the  analytic  satellite  theory  of  Lane  (1965)  and  Lane 
and  Cranford  (.1969)  which  models  the  gravitational  zonal 
harmonics  through  J5  and  models  the  atmosphere  with  a 
spherically  symmetric  power  density  function  is  given  here 
in  its  complete  form.  All  equations  needed  for  satellite 
prediction  are  given  and  the  reader  is  referred  to  the 
original  publications  for  the  theoretical  background.  Two 
simplified  equation  subsets  (IGP4  and  SGP4)  of  AFGP4  are 
also  given  and  the  procedure  by  which  they  were  obtained  is 
outlined.  The  IGP4  and  SGP4  theories  were  originally 
developed  by  Kenneth  H.  Cranford. 
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1. 


INTRODUCTION 


The  analytic  drag  theories  developed  by  the  Office 
of  Astrodynamic  Applications  are  all  descendants  of  a 
parent  theory  AFGP4.  The  lineage  can  be  traced  from  AFGP4 
to  IGP4  and  finally  to  SGP4 .  In  general,  the  theories 
decrease  in  complexity  and  accuracy  in  the  order  AFGP4, 
IGP4,  SGP4.  This  eport  gives  a  comprehensive  description 
of  the  equations  a?  they  are  implemented  and  describes  the 
process  by  which  the  simplified  theories  were  obtained. 

2.  AFGP4 

AFGP4  is  a  general  perturbations  theory  based  on  the 
work  of  Lane  (1965)  and  Lane  and  Cranford  (1969) .  The 
gravitational  model  includes  the  effects  of  the  first 
five  zonal  harmonics  of  the  Earth.  The  atmospheric  model 
assumes  a  static,  non-rotating,  spherically  symmetrical 
atmosphere  whose  density  p  can  be  described  by  the  power 
function 


where  the  quantities  on  the  right  side  are  defined  in 
the  list  of  symbols  in  the  Appendix.  A  discussion  of  the 
development  of  the  equations  is  thoroughly  covered  in  the 
publications  cited  above  and  will  not  be  repeated  here. 
The  numerous  symbols  used  in  the  following  equations  are 
defined  in  the  Appendix. 
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Given  the  values  at  an  epoch  t  of  the  mean  orbital 
elements  (denoted  by  subscript  o)  predictions  with  AFGP4 
are  made  according  to  the  following  scheme: 

a.  Calculate  the  constants 

k  =  By(l-n2)'7/2  {(1  +  |  e2  +  |  n2  +  3eZn2  +  4en  +  en3) 

+  (£)  C-  |  +  1  e2)5(i-n2)-1  [(8  +  24n2  +  3n4) 

~  5en(4  +  3n2)]} 

A,  -  L’V 

A2  =  4  C9  1 

A3  =  (4/3)  L'1 (17a+s) 1 
A4  =  (2/3)  (221a+31s) £19y"2 

where  all  quantities  on  the  right  side  are  understood  to 
be  double-primed  epoch  quantities. 

b.  Calculate  the  time-dependent  quantities 

Q1  [k2  (A2  +  2Al)  Ct-tQ)3  +  0.25k3  (3A3  +  12AxA2 

+  lOAj’)  (t-tQ)4  +  0.20k4  (3A4  +  12A1A3  +  6A2  +  30A2A2 

isaJ)  (t-to)5] 


4* 


Q2  =  -  [A2k2(t-t0)2  +  A3k3(t-to)3  +  A4k4(t-to)4] 

k  k2 

Vi  =  Vo  +  {1  +  I  ~TTS  C‘1  +  3e2)  +  fl  -7TJ  t  ~  1 5  +  166  +  2 56 2 

ci  P  <1  P 

k 

+  (30-966-90B2) 02  +  (10S+144P+2582) 04 ]  +  -i-w  e2  (3-30  02 

10  a  c 7 

+  35  e4)}  n  (t-t01 

k  k2 

*i  -  8o  *  <•  I  TTT  u-se2)  ♦  A  -A  l'55  *  246  -  2Je2 

3.  p  3-3 

k 

+  (90-192g-126B2)e2  +  (385  +  3606  +  4532)64l  +  L-  -J-jr  [21-962 

10  a4p 

+  (-270+12682) 02  +  (385- 189g2) 04]  }  n  (t-t  ) 
k  k2 

hV  =  +  (-3  6  +  |  -fr  [  ( -  5  +  12B+902}  0  +  (-35-368-582)03] 

a  8  a  8 

+  I  TT  e(3-702)  (5-382)}  n  (t-t  ) 

4  a488  ° 

where  all  quantities  on  the  right  side  are  understood  to 
be  double-primed  epoch  quantities. 

c.  From  JlV,  use  the  definitions  of  the  angles  to 
calculate  EV  and  XV. 

i  1 
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Then  calculate 


(%  (Xi 

1  t  =  /  /  (Y.T  1  +  2  ey  7 


— 3  +  — V  JT-2  +  ^1',‘  ^  t — 2  +  — 3 
n3  2n4  2  2nz  n3 


(l-n2)3  +  7~^)  J  +  (1-n2)4  ^ 

n-3  2n  8n 


y  J  1  (Y2T'1  -  f  e2y2y2T"3)dXdt 


o  o 


{(1  -&>  Jt-i  *  71  12  J 

2n  2n 


(i  -  n2)JT.3  ]}(iy  -  r;)2 


dt  -  -  SltUUS -  ,9k2  ,13 


-o-!  oj 


azAl  „sin  I  k,sin  I 


3,0 


32k5 


1  U-n2)T*I  io>t,2  (»I 


o' 


015  *7  in  I  -y 

2BT(T-t-l)e°g  a_  { - 3_,0  -  (sin2I  +  12)  I,  T+2[sin  *i 

9n(l-502)2  8k2  1,T  2  1 


sin  g"  -  (gV 


?  2a2 

„,n  .  sin  I  ra  A3,0 

go5  C0S  go)]  +  —2—  l”7^4" 


( 3 02 - 1) ]  I2jT  +  2  t-  J  (cos  2gV  -  cos  2g^)  -  (gV 


g")  sin  2g|J]  + 


aA,  nsin3I  , 

*-“2  I3  t+2  ^"3  (s:*-n 

24k,  •5>  L  J 


sin  3g^) 


(gV  -  g”)  cos  3g”]  +  4?-1  I4>x+2t-T  (coS  4gi  "  C0S  4go} 


[g£  -  gJJ)  sin  4g”]  } 


f  t 


AL . dt  = 

l 


-  Baa  [1  -  (§)  3  *  f  e2)]  Ij 


bx25t  t-l  *  f  e2}  &[1  •  (t)3  (4*1  02„  (tS 


r)Z  {-  t£[  -3 IQ  jT+1  -  I  e  Ij_  >T+1 


+  eZ  ^‘2  *o  ,t+l  ■  1 2  , t  + 1  ^ 


e3  ( -  —  i 
e  1  16  1,T+1 


+  2  (gV  -  g")  sin  2g^] 


5  8 


a  rlrTA,  sin  I  p 

t  5  3,o  j.tM-t  +®r-T 

+  — J—  *2.2  ~  1  TU  2  l  Ao,x+] 


Qk^ntl- 59  ) 


'  3I2,t  +  P  +  I  t_3Il,T+l  ‘  9I3,T+l) 


16  ^3Io,t  +  1  "  8I2,t  +  1  "  11  1 4 ,  t  +  1^  ^ 


*  !  [‘l,t  +  !  00jT  -  3I2,t)  *  b  (9I1,t 


+  11  I^T)  ]  }  [sin  gV  -  sin  g” 


(g  V 


gp  COS  g”  ]  +  6 


f\ 

'  t 


dt  +  Qx 


/  AG.  dt  =  /  AH. 

<.  1  >L 


dt  =  -  Baa[l 


3k. 


+  ?  02^  lZ  ‘  In 


3BkzC 


C-  I  +  I  ^^O.T  +  l 


.  4  ,-T  „  8 


-  J l 


»)2  -  B*X&  I1±4I«  (cos  2g” 

o  18nk2  (1-56Z)2 


-  cos  2g^  +  2 (gV  -  g”)  sin  2g^  } 


1 


8 


§ 


^2  /  t 

+  3n  7171  ( 2  -  158  2 )  )  ag 

t 

o 


a'B 


dt 


2  2 

+  15n  -J-r  d 

a  V  't 


/‘-I 


dt 


hi  +  l 


hi  +  9n  T7 
a  0 


+  1 5n 


k2e 

~ 2~4 
a  a 


o 


dt 


dt  *•  3n 


k26 

a234 


/tiHi 


dt 


where  all  quantities  on  the  right,  side  are  double- 
primed  epoch  quantities  unless  otherwise  noted.  Repeat 
step  c,  using  *V+1,  g^+1>  hV+1  in  place  of  AV,  gV,  hV , 
respectively.  The  number  of  iterations  required  is 
determined  by  the  value  of  X  =  |  n  dt  | 

according  to  the  following  table: 


X 

Iterations 

[ 

!  ,5  <  X 

1  — 

4 

j  ,05  <  X  <  ,5 

2 

* 

A 

O 

1 

e.  Let  i's ,  g’s',  hg  denote  the  final  iterated  values, 

and  let  A”  and  E"  denote  the  values  of  A"  and  E"  corre- 
s  s 


sponding  to  V.  Then  calculate 


I 

J 


m  sinX  ( 2  y2  T  ~ 2  +  2eYlY2T'3NX  =  -  ^  Id 

2 

+  n  cos  X^)T  1  -  (1  +  n  cos  X")T  *]  •  l  (1 

+  n  cos  X'pT  2  -  (1  +  n  cos  X^)T  2] } 


•  T-l, 


Yttt  l  Cl  +  n  cos  \'p  '  -  (l  +  n  cos  a"j  ] 


1  P  BT(Tnkl^  [9k2  (.  1  +  3  92)  2  +  sin2l 

L  2  2a  III  32k^ 


->  4 

k.,sin  I 

♦  ^-T~  ) 


fX%2^idx  ,  BT.^.niV^g. 

)  ,  2  3  (1  -  50  2  ) 


aA,  nsin  I  7 

•-  — -  f  sin2  I  +  12)  1 


(sin2  I  +  12)  l1>T+2tC0S  £5  "  cos  £o] 


.  2  T  a  A  *7  n  3k r.  ~ 

+  MS_I  [ - ^0  +  — £  t3e  2  -  1)  ]  I 

2  32k^  4 


2  ,t  +2  [sin  2gs  '  Sin  Zgo  1 


10 


■V  ! 


q  ^4  -?  oft4 

'  36  JZ  *  1'36  '  (sin  2Ss  -  Z$") 


A 

J  “7i'a  sin  1  (cos  gs  ■  cos  S”)>  t  -  TCI  -  I-L  . 

J\  o  * 


2  e  ^Ojt+l  +  3I2,t  +  i^  +  a  i  r’2^  +  Lq  ^2 


[  1<-Baa[l 


10 


(•-7*1  »2)) 


7  e2i!,r2r'3) dx  -  ek2{*  (.1  .  |  92)  U" 


-  tp  I  re  I-  31^^  *  |  eli>T,j]  *  \  13I0it!) 


+  Ba2gTB4  (l-e2-) 
1-5  0  2 


{-tC  [  jr  I 


-  I  n 


6  2  ,t  +  1  "  6  1  1,T  +  1  r  3,xtl 


+  Ji  I2,T>[sin  2«S  '  sin  2*o] 


! 


♦  (  i  [  1  - 1 ie 2  -  foe_] 

3(1-50  8  1-50* 


[1-  302  ‘  “iTflHxC  [eICjT  +  1  ♦  eI2>T  +  1] 


a  I1,t}  tsin  2gs  '  sin  2go] 


BB6a3£TA,  sin  I 

- ^7,-2-  *  '2,,. 


*  e  h.fl11  *  5  (10,T*l))l<:05  8S 


-  cos  g”] 


and  the  drag  periodics 


e”5*D  =  Baa Y  1-n2  I3  +  BaVfi5  ^  {-  \  ulI3>T  +  1 


-  1 1 ,  t  + 1  +  e  c  i 1 


1  1  _  25  T  +  15 

4  o,T+l  T  ^,x  +  l  T  4,t+1 


l  h'T  +  l»  +  h  l-3Il,T  *  7I3,x1  }  tcos  2gs 


COS  7,g"  ]  +  {  2  Cl- lie2  -  ~^~~T  ) 
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'D 

where  all  quantities  on  the  right  side  (except  g” , 
E'g'  and  X")  are  understood  to  be  double-primed  epoch 
quantities  and  where 


Jy^y^dy  YL  TlP"k+S'^)  (J)  cosk+WAsinA 


k  =  0  4-0 


i  - 1 

+  H  . 

j  =2  l 


k=2  (^k  +  3)  __  k+4-2j  +  1\  _  . 


cos 


n 

k=l 


(4-k+2) 


A  s  i  n  A  +  R . } . 

l  ’ 


k+4+1 


if  k+4  is  odd 


with  i  H 


k+4+2  . r  .  o  . 

—  if  k+4  is  even 


and 


^A  if  k*4=0 

0  if  k+4=l 

I  7  k  -  A 

II  ( 2 k -~2 ^  ^  k*-'  *s  even  ar>d  k+^2 

k — 2 

1  7  v  _  7 

h sinA  if  k+4  is  odd  and  k  +  &>3 
k=2ZK_i 
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V. 


a"  “  tan  1  [e"  sin  s."  /  e"  cos  t"] 

h"  =  tan'1  [(sin  I**/2  sin  h“)  j  (sin  l"/2  cos  h”j] 

g"  =  U"  +  g"  +  h")  -  £"  -  h” 


I"  =  2 
and  L" 


sin'1  V(si^ 
has  already 


I"/2  sin  h")2  +  (sin  I"/2  cos  h" 
been  given  in  step  e. 


g.  The  periodics  due  to  the  geopotential  model 
now  included  by  first  calculating 


n  *  U"  +  g"  +  h")  + 


4a2g2  Cl^)  U('1  +  3e2) 


*  f  +  1)  sin  f  +  3  ( 1  -  0 2 1  [(-  .  a  + 

r2  r 

+  11  sin  C2g  +  f)  +  (^f2  +  f-  +  1 )  sin  ( 2g 
3  ^2  ? 

+  2  ~77?  (-  1-29  +56^)  (f  -  £  +  e  sin  f) 
a  p 

k2 

+  4a"5 g4  C3  +  2 0- 5 02D  [3  sin  (2g  +  2f) 


+  3e  sin  (2g  +  f)  +  e  sin  (2g  +  ^ ) ] 

+  IE  ~T7$  {(203-2-e2)  -  11  (233-2-3e2)02 
a  b 

18 


+  3f)]} 


k _ t 


vtfrm  "I  ffHW‘111  WW’HU  K»ii<lHMi».BIl!l|[.|ff lilpiyWII^llM^^ 


-  4O(283-2-5e2)04  (i-502)'1  -  400  e206  (1-502)"2 


-  2e26 [ 11  +  8062  (1-S62)'1 


+  2000 4  (1-50  2)  " 2  ] }  sin  2g  -  Jl  -  -A— a  {(2B5 


V  B 


-2-e2)  -  3  (2B3-2-3e2)02  -  8  (203-2-5e2)  — - 7 

1-56 


+  8Oe206  (1-56 2)  2  - 2e2  0 [ 3  +  1602  (1-502)*1 


A 

+  400 4  (1-562)"2]}  sin  2g  +  i  544  [  (1+B 

V6 


+  02)  (— ^— )  sin  I  +  e9  tan  1/2]  cos  g  + 

1  +  B 


+ - .  [402  Sin  i  (  ®  )  +  3eB2  sin  I  -  9eB3  sin  I 

k2a-V  P 


+  (4  +  3e  je  0tan  1/2  +  e  sin  I  (26 


+  9e2) ]  [ 1  -  90 2  -  2494  (1-502)"1]  cos  g 


+  44  5£- ^  e0  sin  I  (1-0)  (4  +3e2)  [3  +  1662  (1 


-56  )  +  408  (1-50  )  ]  cos  g 


+  -^Lr  - |2_  [e  sin  I  (3B3  -3  - 2e2) 


1152 


•e3  Gtan  1/2]  [  1  -  56 2  -  16B4  (1-592)'1]  cos  3g 


35  50  3  .  2  2 

5^  - fL-  e-5  Asm  I  (1-0)  [ 5  +  326 ^  (1-50^) 


k2a  8 


+  8O04  (1-502)  2  ]  cos  3g 


L"  U  +  — [3(f)3  (1-02)  cos  ( 2 g  +  2f)  +  (362 
2a“ 


i)  ((f)3  -  e~3)]} 


4a2B 


2  2 

2  2(-l  +  302)  (^-  +  f  +  1)  sin  f 


?  2r2 

+  3(lV)  [(-  -  f  +  1)  sin  ( 2 g  +  f) 


2  2 

+(Mr  +  I +  sin  (2g  +  3f)1} 


+  {^C  f  1  -  1 16 2  -  4 0 0 4  (1-502)'1] 


9 


^4e  ,,  ,„2  ori4  n.ca^'ll 


1  ^3  O8 

+  -  4-  4^-  sin  I 


_ !_  [1-30^  -  80  ^  (1-56  )'  ]}  sin  2g  +  -  4  “kT 

12k2a  ^ 


.  S  sin  1  (4  +  9e2)  [  1  -  962  *  2404  (1 


64  k2a303 


Se2)"1])  cos  s  *  OT  — ^#3  e2  1  fl-S8Z-l«e2  (1 

384  k2a-iB'5 


- 50 2 )  1]  cos  3g 


fie  =  fi ,  e  +  — 2-4  { (-l  +  302)  [e&  +  yfg  +  3  cos  f 
2  a  8 


+  3e  cos^  f 


+  e2  cos3  £]  +  3(1-4)  (e  +  3  cos  f  +  3e  cos  f 


+  e2  cos3  f]  cos  (2g  •  2f ) } - 4y  (1-62)  [3  cos  (2g 

26  a 


+f)  +  cos  (2g  +  3f)] 


where  626  =  sin  I  fi?e»  where 


6  e  =  {1  44  sin  -  —-*-TT2  {ll23~l  Sin  1 

2  8  a202  1  -  56 2  12k2a  6  1-50 


]}  cos  2g 


+  {I  4r°+  JL  ULiO  (4  +  3e2)  [l-902-2404  (l-502)  1]}  sin  g 


4  V  64  k2aV 


f  A5,0  - 

(vV 


JJL  ^  As  .Q  )  e2  [  1- 59 2  - 1604  (1-502)'1]  sin  3g 


i.  l)i<  Klliivij  JJ  il  Ilk!  ilijuikL  >  III  li  kiiii  i  ^ :  lllUjlliliL:!  :iil  )il  iillililikllliikl  liiMiji  I  llLlliliy  ti  nil  li 


+  44  - tTT  e  6sin  1/2  sin  I  (4  +  3e  )  [3 

3  L  i  _  3  0  o 


V  e 


+  i6e 2  (l-se2)’1  +  40e4  (i-se2)'2]}  cos  g 


35  A50 


e3e 


+  t-  ifti  7^X6  (  r^iTTT)  t1'59 


V  * 


i6e4  (i-se2)"1] 


- -3P,  e3  esin  I  sin  1/2  [5  +  3202  (1 


576  k2a306 


-502)  4  +  8O04  (1-S02)  2]}  cos  3g 


H  =  H" 

where  all  quantities  on  che  right  side  are  understood  to 
be  double -primed  elements  at  time  t. 

h.  The  osculating  orbital  elements  can  now  be 
determined  by  first  calculating 

sin  1/2  cos  h  =  cos  h"  [(61/2)  cos  I"/2  +  sin  I"/2] 


sin  h"  sin  (I"/2)  <$h 


cos  1/2  sin  h  =  sin  h"  [61/2  cos  I"/2  +  sin  I"/2] 

+  cos  h"  sin  (I"/2)  6h 
e  cos  i  =  (e"  +  6e)  cos  £"  -  e"  6£  sin  £" 

e  sin  SL  =  (e"  +  6e)  sin  £"  +  e"  6SL  cos  X," 

i  +  g  +  h  =  £"  +  g"  +  h" 

Then,  the  complete  set  of  osculating  elements  at  time  t  is 
e  =  V  (e  sin  +  (e  cos  £)2 
&  =  tan  1  fe  sin  £/  e  cos  2-] 

h  *  tan  1  [(sin  1/2  sin  h)  /  (sin  1/2  cos  h) ] 

g=(£+g+h)-£-h 

1=2  sin  1  V  (sin  1/2  sin  h)2  +  (sin  1/2  cos  h)2 
a  =  L2/p 
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3 .  Drag  Simplification 

The  drag  portion  of  AFGP4  is  simplified  by  using  only 
the  secular  terms  from  the  equations  along  with  certain  drag 
periodics  of  important  magnitude.  The  value  of  t  is  fixed 
at  4  throughout  the  simplified  equations,  and  all  quantities 
are  understood  to  be  mean  epoch  elements  unless  otherwise  noted, 
a.  The  L"  equation  from  AFGP4  (section  2e.)  is 


A"  •  A” 
s  o 


*s  -  £o  -  no  (t  •  V 


and  drop  certain  smaller  terms,  we  have 


L"=  L”  [1  -  Cl  (t  -  tQ)]  +  L”  Q2 
where 

C1  =  B£4n  iJ/~7  [a  (1  +  |n2  +  4en  +  en3)  +  \  <J>'2  (-  \ 

+  |92)  (8  +  24n2  +  3n4)] 

Dropping  the  same  small  terms,  the  constant  k  (section 
2a.)  is  approximately  by 

k  =  na2C^£  4 

It  then  follows  that  Q2  (section  2b.)  can  be  written 

Q,  =  -  D,  (t  -  t  )2  -  Dt  (t  -  t  )3  -  D.  (t  -  t  )4 
x2  2  v  o  3  o'  4  v  o' 


where 

Dz  -  4a  £  Cx2 

U3  =  Ja^2  (17a  +  s)  Cl3 
d4  =  |  a2?3  (221a  +  31s)  C ^ 


Then 
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L"  =  L»  [1  -  Cx  (t  -  tQ)  -  D2  (t  -  tQ)2  -  D3  (t  -  to)3 


'  °4  (t  '  to)  5 


b.  From  section  2c.  ,  we  make  the  approximation 


/ 


ALi  dt  =  -  Baa  ^  -  Bk^4  (-  \  +  jQ2)  U- 


-  ft.”)2  12U0)S  +  Qj 

Using  the  same  approximation  as  in  section  3a.,  we  take 


ALi  dt  = 


\  C1  (t  -  t0)2  +  Q1 


t 

o 

where 


Ql  =  '  |  t (°2  +  2  Cx2)  (t  -  tQ)3  +  \  C3D3  +  12  Cx  D2 

+  10  C13)  (t  -  tG)4  +  |  (3D4  +  12  C1  D3  +  6D22 

+  30  Cj2  D2  +  15  Cx4)  (t  -  to)S] 

It  is  further  assumed  that 
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..'/jl  |Mf 


The  integral  I-  is  approximated  by 


h  -  -  H  f11*11  cos  ii) 


»•  ’i  ^  _ 


(l-n  cos  ^)3] 


so  that 


6*„-  -  |  BacS'6  CCS  tv)3  -  (l-n  cos  t”)*J 


+  BaBn£Sifr"9  j^0  sin  I  (l  +  3nZ  +  |n4)  cos  g  tt 


where  the  approximation 


sin  g»  •  sin  g”  -  (t  -  tQ)  ^  (sin  g’’) 

—  k  ^  n  j 

=  -  4  -4s  u-50  )  sin  b  (t  •  t0) 

1  a 

has  been  made.  Finally  with  further  approximations, 
we  take 

6Hr  «  -  I  BC4  — r  [  (l+n  cos  Z'P5  -  (1+n  cos  fc”)3] 


+  B*  C3  cos  g  (t  -  tQ) 


where 


S  ■  <% 


s)4{5  an  3^  Sin  I 
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I 


The  drag  period  Ag,,  is  taken  to  be 


%  "  -«*D 


d.  The  e”  equation  from  AFGP4  {section  2e.) 
approximated  by 


eD  =  eo  ~  Baa®' 


(2YiY2  +  2eY1Y':Y2)dX 


BS255k2  (■  1  *  |s2)  (E”  -  E^)  | 24  Ilj5 


’  6e  (Io,5  *  '2,8)  *  if  ‘1,4) 


♦  2  BSV55  [■  j  I,  ,  •  I  I,  , 

c 1 -  5  e z )  3  1'5  3  3>s 


+  e  (4  I  c  +  4  I 


■3  o,5  .5 


I  X 

Z , 5  '  3  ?5^  +  2a£ 


+  I  h,V]  (sin  2gs  •  sin  2go} 


The  integral  term,  when  evaluated,  is 


-  2  B^’7?4  aB2  (2n  +  In3  +  \  c  +  2en2)  (A"  -  A") 

z  z  o 


-  2  “ 7 C 4 a B 2  [(1  +  +  ^en  +  en3)  (sin  X" 


sin  X^)  +  (yn  +  jn3  +  ^-e  +  (sin  2X' 


-  sin  2X^)  +  (^-n2  +  j^en)  (sin  3X"  -  sin  3X^)] 


which  is  approximated  by 


2  B^'7£4a82n  (2n  +  |n3  +  \e  +  2ep2)  (t  -  tQ) 


■  2  B^‘7£4a$2  (1  +  Q n2  +  ^en  +  en3)  (sin  A' 


-  sin  A”) 


where  1"  =  A”  +  <5AD.  The  second  term  gives 


4  B62C5ip‘9k2n  ( 1  -  30  2 )  (  ^n  +  ^-n3  +  §«  +  (t  -  tQ) 


where  the  approximation 
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E 


E"  - 

S  0  O'  O' 

has  been  made. 

The  third  term  given 

2B32cV9k2n  Cl- e2)  (ion  +  ^n3  -  J  e  -  13en2 

+  cos  2g 

where  the  approximation 

sin  2g^  -  sin  2g^  =  (t  -  tQj  ^  (sin  2gs"] 

=  -  3  (1-502)  cos  2g  (t  -  tQ) 

a  E 

has  been  made.  Thus,  the  final  equation  for  ejj  becomes 

eii  -  e0  '  B*  C4  -  to)  -  B*  C5  lsin  <*s  +  5V 

-  sin  it") 
o 


wheret 

C4  =  2  (qo  -  s)AZ  V7  a  B2n  f  [  2  n  +  |n3  +  Ie  +  2en2] 


t  It  should  be  noted  that  the  k2  terms  of  in  the  opera 
tional  version  do  not  agree  with  this  equation.  The  pres 
ent  equation  has  been  verified  as  correct  with  the  author 
and  theoperational  version  will  be  modified  to  agree. 
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-  2  ^  r2[(l-3  02)  ( 

-  ^en4)  +  (l-62)  (Sn 
+  ^eb4)  cos  2g]  > 


IS  J5  3  +  9  9  2 

Tn  +  Xn  Ie  4*n 


.  65  3  7  13  2 

+  -rn  '  ^ 


e . 


(qo  -  s]4  €4i|>"7a62  (1  +  -^n2  +  ^en  +  en3) 
The  double-prime  updated  elements  are 


I"  =  I" 
o 

*"  =  *s  +  6% 

g”  -  g5+  sgD 

h"  -  h" 
s 

and  L"  has  been  given  in  section  3a. 
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4. 


IGP4 


IGP4  is  obtained  from  AFGP4  by  using  the  drag 
simplification  discussed  in  section  3  while  retaining 
the  full  Brouwer  geopotential  (1959)  with  the  Lyddane 
jnodif ication  (1963)  . 

Given  the  values  at  an  epoch  t  of  the  mean  orbital 
elements  (denoted  by  subscript  o)  predictions  with  IGP4 
are  made  according  to  the  following  scheme: 

a.  Calculate  the  constants 

C2  1%  ~  s)  5  n  u  "  h  )  [a  (1  -  n  )  (1  +  2-n 
4en  +  en3)  +  §  (-  \  +  ffl2)  (8  +  24n2  +  3n4)] 


B*C- 


C,  = 


(qo  '  s)4C5A3  Qn  a  sin  I 
k2e 


c4  =  2n  (%  '  s^>454aB2  (l  -  nV7/2([2n  (1  +  en) 


(■ 


13, 


+  2e  +  2n 


2k2? 


a  (l  -  n  ) 


t  f,  7  a 2 s  ,  15  .  45  3 
2  [  (1-36  )  (  -y-n  +  -g-n 


9  9  2  3  4s  A  ,,  „2s  .  65  3  7 

+  2i  4en  "  T )  +  U-e  )  (5n  +  -5-n  -  ^ 


IF 
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-  iren2  +  ^eq4)  cos  2g^j 

c5  =  2  (qo  -  s)Vae2  (1  -  n2) " 7/2  11  +  Yn  (n  +  e) 

+  en3] 

D2  =  4a  ^C12 

D3  =  |  a  (17a  +  s) 

D4  =  |  a  2C3  (221a  +  31s)  q4 


where  all  quantities  on  the  right  side  are  understood 
to  be  double-primed  epoch  quantities. 

b.  Calculate  the  time  dependent  quantities 


r.' 

l 


*o  *  U  *  !  -T7s  (-1+3eZ)  *  A 


a  3 


15  +  163 


+  25g2  +  (30  -  96g  -  9032)B2  +  (105  +  1443 


+  2S32)94]  +  ~  — 4~ y  C 3 - 3 0 6 2  +  3 5 6 4 ) }  n(t  -  tQ) 

a  3 


go  +  {'  I  -!r?  C1'502)  +  ji  4  t"  35  +  24B 

O  A  J>- 
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+  25g2  +  (90  -  1923  -  126g2)62  +  (385  +  3603 


+  4Sg2)04]  ♦  yk  -Ag  121  ‘  9b2  +  270  +  126B2)92 

a  3 


+  (385  -  18932)94]}  n(t  -  t  ) 


k  k2 

hi  "  ho  *  <  -  3  TI74  6  *  1  T78  l(-  5  *  126  ’  9s2)e 


a  3 


a  3 


+  (-  35  -  363  -  S32)93]  +  |  -J-~  9  (3-792)  (5 

a  3 


-  33  ) }  n (t  -  tQ) 


where  all  quantities  on  the  Tight  side  are  understood  to 

be  double-primed  epoch  quantities. 

c.  Let  Si'',  g" ,  h"  denote  the  secular  values  of  l", 
s  7  °s  s 

g" ,  h" ,  respectively.  Then 


g”  “  g 


h" 

s 


h'.' 

i 


21  nk26 


2,2  °1 


a  3 


(t 
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*s  =  ri  *  ”  if  C1  <l  •  ‘o’2  *  (°2  *  2C!2)  (t  -  t0)3 

-  i  (3D3  *  12CxD2  *  lOCj3)  (t  -  to)4  -  |  (3D4 

+  12C1D3  +  6D22  +  30C12D2  +  lSCj4)  (t  -  tQ)5] 

Also 

L"  •  Lo  I1  •  C1  Ct  -  t0)  -  D2  (t  -  tD)2  -  D3  (t  -  t0)3 

-  »4  (t  -  t0)4l 

eD  =  e0  *  b*C4  (t  -  tQ)  -  B*C5  [sin  (*”  +  6Ad) 

-  sin  l"] 

I"  -  I" 

0  o 

where  all  quantities  on  the  right  side  are  understood 
to  be  double-primed  epoch  quantities, 
d.  Calculate  the  drag  periodics 

“d  *  -  §  <%  -  S^B*54  £  Hi  *  n  C05  tv)3 

-  (1  +  n  cos  A")3]  +  B*C,  (cos  g")  (t  -  t  ) 


«8d  -  -  .&lD 

where  all  quantities  on  the  right  side  (except  &V) 
are  understood  to  be  double -primed  epoch  quantities. 

e.  The  complete  set  of  double -primed  elements  at 
time  t  is 


e 


11 


i"  =  j" 

D 


%"  = 


6£ 


g"  -  g^  +  68D 


and  L"  has  been  given  in  section  4c. 

f.  The  geopotential  transformation  for  IGP4 
osculating  orbital  elements  is  then  given  by  the  equa¬ 
tions  of  section  2g  and  2h. 


3  8 


. . »- . I'fBiupi'iwi  . . . HI .  n»i  . . . . . . --nr 
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5.  Geopotential  Simplification 

The  geopotential  portion  of  AFGP4  is  simplified 
by  assuming  eccentricity  is  small  and  then  retaining 
only  the  main  terms.  The  small  eccentricity  problem 
is  handled  by  using  an  alternate  set  of  variables. 

The  small  inclination  problem  is  not  treated  since 
few  satellites  which  experience  drag  forces  have 
snail  inclinations. 

a.  The  secular  geopotential  equations  (section  2b.) 
are  simplified  by  taking  £-1  except  for  those  £  appear¬ 
ing  in  the  denominators.  The  result  is 


q  =  q  +  [i  + 


3k,(l-  +  36') 

_7? 


3k22  (13  -  7 8 8 2  +  1 3 7 0 4 


16a43 ' 


i  n(t  -  t 


S'l  =  Sq  +  I- 


3k,  (1  -  se2)  3k,2  (7  -  114e2  +  39504) 


tv  -  360“  +  490') 


5k40  (3  -  78“) 


]  n(t  -  tQ) 


where  all  quantities  on  the  right  side  are  understood 
to  be  double -primed  epoch  quantities. 

b.  The  long-period  terms  of  the  geopotential 
(section  2g.)  are  approximated  by  retaining  only  the 
main  Q  term.  The  result  is 


6(1"  +  g"  +  h").  +  4 


1  A3,0 
8  k2a$2 


lin  1  tjTe  )  cos  R 


{Ll  ’  0 


«L  '  0 


*hL  '  0 


6e.  =  i  -  3 ’  —j  sin  I  sin  g 

1  4  k2ar 


13  0  1 

6h  =  4  e  Sin  1  cos  g 

k2ae 


where  all  quantities  on  the  right  side  are  understood 
to  be  double-primed  elements  at  time  t,  and  where  the 
subscript  L  denotes  long-period  geopotential  periodics 


V 

t 

I ' 


\ 


>i  nib 


If  we  consider  the  two  new  variablese’  sin  g’  and  e1  cos 
g’  ,  then 

- e1  sin  g1  =  e"  sin  g"  +  6(e"  sin  g")^ 

=  e"  sin  g"  +  6eL  sin  g"  +  e"  fig^cos  g" 

e'  cos  g'  =  e"  cos  g"  +  6(e"  cos  g")^ 

=  e"  cos  g"  +  6eL  cos  g"  -  e"  <5gL  sin  g" 
Upon  substituting,  we  find 

1  0 

«5(e"  sin  g")  =  4  -  sin  I 

kza6Z 

6(e"  cos  g")  =  0 

Then  the  long-period  geopotential  transformation  is 
V  +  g*  +  h'  =  %"  +  g"  +  h"  +  6  (&"  +  g"  +  h") L 

e'  sin  g'  =  e"  sin  g"  +  S(e"  sin  g")L 
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fc,  i 


e'  cos  g'  =  c"  cos  g" 

L'  =  L" 

I*  =  i" 

h’  =  h" 

and  it  follows  that 

e'  =1 1  (e '  sin  g' )2  +  (e'  cos  g * ) 2 

c.  Before  adding  the  short-period  terms,  a  change 
of  variables  is  made  according  to  the  following  scheme. 
Let  E'  denote  the  eccentric  anomaly  corresponding  to  . 
Solve  Kepler's  equation  for  E'  +  g'  (by  iteration  to  the 
desired  accuracy)  where 

(E*  +  g’)i  +  1  =  (E'  +  g')i  +  A (E '  +  g')i 


with 


A(E '  +  g')i  = 

£'+g'  -  (e'  sin  g')  cos  (E'+g')^+(e'  cos  g')  sin  (E'-t-g')^  -  (E'+g'^ 
-  (e*  sin  g')  sin  (E'  +  g')^  -  (e'  cos  g')  cos  (E'+g')^+l 
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and 


(E '  +  g')x  =  +  g’ 

Then,  compute  the  following 

e*  cos  E*  =  (e‘  cos  g')  cos  (E '  +  g’) 

+  (e '  sin  g' )  sin  (E’  +  g’ ) 

e’  sin  E '  =  (e‘  cos  g’)  sin  (E'  +  g’) 

-  (e1  sin  g')  cos  (E ’  +  g’) 

p'  =  a'  (1  -  e' 2) 

,  ,  _  1  , ,2 
where  a  -  —  E 

and 

r'  =  a'  (1  -  e'  cos  E’) 

ri  =  kl  e'  sin  E' 
r ' 

,  •,  .  dJ-EL 

where  f’  denotes  the  true  anomaly  corresponding 


to 
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Let  u'  =  f '  +  g' .  Then 


n  * 

cos  u'  =  [cos  (E '  +  g'J  -  e'  cos  g' 


+  (e1  sin  g1)  (e1  sin  E1)  ^ 


1  +' 

rr 

,  2 
-  e' 

sin  u 1  —  —y 
r  * 

[sin 

(E* 

+  g' )  -  e' 

(e1 

cos 

JC2 

(e '  sin  E ’ ) 

1  +1  1-e'2 


At  this  point  we  have  the  following  set  of  single-primed 
variables . 

•  • 

r' ,  r' ,  r ' f ' ,  I' ,  h' ,  u’ 

d.  For  the  above  set  of  single-primed  variables,  we 
can  write  the  following  formulas. 

(<5r)s  =  ^  sin  f  (6Z)S  +  |  (6a)g  ’  a  cos  f  (<$e) s 

(6r)s  =  -  \  ^  sin  f  ( Sa)  s  +  ~  (|0  (1  +  e  cos  f)  sin  f  (<Se) 

+  nae  cos  f  (^) 2  (6£)s 
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I 


/  V 

t  '  j 

I .  J 

■  j 


6(rf)s  =  -  i  ne  (|)  (6a) s  +  ^  (|)  (cos  £  -  e 


+  e  cos2  £)  ( 6e)  -  nae  sin  f  (^)2  (6H) 


(6u)  =  -L  (2+  e  cos  f)  sin  f  C<Se)  +  -V  (1 

s  e2  6 


+  e  cos  f)  (6*)s  +  (6g)s 


where  all  quantities  on  the  right  side  are  single-primed 
quantities  and  when  the  subscript  s  denotes  short -period 
geopotential  periodics.  We  then  substitute  from  the 
Brouwer  formulas,  simplify,  and  make  the  approximations 


(6r)s  =  \  (1-e2)  cos  2u 

a  (3 


(6r)c  =  0 


6 (rf ) s  -  0 


-  k? 

61  =  j  ■  2  '4'  6  sin  ^  cos 

a  g 


1; 

6h  =  j  ■■  6  sin  2u 

a"e 


I  (302- 1)  r 
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6u  =  -  i  — ^  ( 76 2 - 1 )  sin  2u 

where  all  quantities  on  the  right  side  are  single-primed 
quantities. 

e.  Osculating  variables  are  given  by 

r  =  r’  +  (Sr) 

s 

•  • 

r  =  r  ’ 

•  • 

rf  =  r ' f ' 

I  =  I'  +  (<5I)S 

h  =  h’  *  (Sh)s 
u  =  u'  +  (6u) s 

and  provide  a  very  convenient  set  of  variables  for  calcu 
lating  position  and  velocity  in  Cartesian  coordinates. 
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6.  SGF4 


SGP4  is  obtained  from  IGP4  by  using  the  geopotential 
simplification  discussed  in  section  5. 

Given  the  values  at  an  epoch  t  of  the  mean  orbital 
elements  (denoted  by  subscript  o)  predictions  with  SGP4 
are  made  occording  to  the  following  scheme: 

a.  The  drag  and  secular  portions  of  the  equations 

are  given  by  following  the  equations  of  section  4a  through 
4e  with  the  single  exception  that  the  secular  geopotential 
effect  of  section  4b  should  be  replaced  by  the  shortened 
secular  geopotentiai  effect  given  in  section  5a. 

b.  The  long-period  geopotential  periodics  are  then 
added  to  the  double-primed  elements  at  time  t  in  the 
following  manner. 

Compute  the  long-period  gravitational  periodics 


6U"  +  g"  +  h") 


A3,0  s1” 
8  k  2  a 


I 


,  ■,  r  3  +  50 , 

(e  cos  g)  (x  T  6-) 


A,  n  sin  I 

5(e"  sin  g")  =  - * - 2 

4k2ag 

where  all  quantities  on  the  right  side  are  understood 
to  be  double -primed  elements  at  time  t. 

Then 


g'  +  h' 


+  gM  +  hM  +  +  g"  +  h”) 
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e'  sin  g'  =  e"  sin  g"  +  6 (e"  sin  g") 

e'  cos  g'  =  e"  cos  g" 
it  follows  that 

e '  ="V  (e~  sin  g ' )  ^  +  (e^  cos  g '  )  ^ 

Let  E'  denote  the  eccentric  anomaly  corresponding  to  a ' . 

Solve  Kepler's  equation  for  E*  +  g'  (by  iteration  to 
the  desired  accuracy)  where 

(E'  +  g')i+1  =  (E*  +  g’).  +  A(E'  +  g'). 

with 

A  (E  '  +  g')i  = 

A'+g'  -  (e'  sin  g')  cos  (E'+g')i  +  (e'  cos  g')  sin  (E'+g')i  -  (E'+g')^ 

-  (e *  sin  g')  sin  (E'+g')^  -  (e'  cos  g')  cos  (E'+g')^  +  1 

and 

Ce  '  +  g'Dj  =  a'  +  g' 

Then,  compute  the  following 

e'  cos  E'  -  (e'  cos  g')  cos  (E'  +  g')  +  (e'  sin  g')  sin  (E' 
+  g') 
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where  £'  denotes  the  true  anomaly  corresponding  to 
Let  u'  =  f ' +g' .  Then 


iniqjlIlia^iiiilliMilMilfiJtULillUtllllliyi^  fcjliiliial  dillillial  iBhik..  Julkilla.. 


u'  =  fr  [sin  ( E '  +  g  ’ )  -  e’  sin  g’ 


(e '  cos  g '  (e 1  sin  E1)  ^ 


1  +  r  1  *  e' 


At  this  point  we  have  the  following  set  of  single-primed 
variables 


r* ,  r' ,  r'f ' ,  I ' ,  h* ,  u' 


c.  Compute  the  short-period  gravitational  periodics 


7  7  \  ^7  7 

6r  =  25  (1'e  J  cos  ?-u  *  j  -fy  (38  -l)r 

p  a  3^ 


6u  = 


\  ( 7 6 2  - 1 )  sin  2u 


3k?8 

6h  =  — =—  sin  2u 
2P 


3k29 

61  =  — sin  I  cos  2u 

2p 


where  all  quantities  on  the  right  side  are  understood  to 
be  single-primed  variables  at  time  t. 


f 


l 


d.  The  osculating  variables  are  now  given  by 

r  =  r'  +  Sr 

u  =  u'  +  6u 

h  =  h’  +  6h 

I  =  I’  +  61 

«  « 

r  =  r' 

•  • 

rf  =  r'f1 

where  all  quantities  on  the  right  side  are  understood  to 
be  single-primed  variables. 
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Appendix 


a  =  semimajor  axis 
e  *  eccentricity 
I  =  inclination 

y  =  product  of  Newton’s  gravitational  constant  and  the 
mass  of  the  Earth 
L  --y  ya 

G  =  L  V~l-e? 

H  =  G  cos  I 

£  =  mean  anomaly 

g  =  argument  of  perigee 

h  =  longitude  of  ascending  node 

a"  =  "mean"  semimajor  axis 

e"  =  "mean"  eccentricity 

I"  =  "mean"  inclination 

L"  =-nri* 

G"  =  L"  1/  1  -  e"2 
H"  =  G"  cos  I" 
l"  =  "mean"  mean  anomaly 
g"  =  "mean"  argument  of  perigee 
h"  =  "mean"  longitude  of  ascending  node 
=  secular  value  of  s, " 
gg  =  secular  value  of  g" 
hy  =  secular  value  of  h" 
qo  =  geocentric  reference  altitude 
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pQ  =  atmospheric  density  at  q0 

s  =  fitting  parameter  in  density  representation 
r  =  fitting  parameter  in  density  representation 
(restricted  to  integral  values) 

E"  =  "mean"  eccentric  anomaly  and  has  the  same  func¬ 
tional  relation  to  e"  and  £"  as  the  eccentric  anomaly 
E  has  to  e  and  % 

f"  =  "mean"  true  anomaly  and  has  the  same  functional 
relation  to  e"  and  2,"  as  the  true  anomaly  f  has 
to  e  and  £ 


n" 

11 

e" 

-  s 

cos 

A" 

cos  E"  -  n" 
1-n"  cos  E" 

sin 

A" 

J  l-n,,Z  sin  E" 
1  - rj ' ’  cos  E" 

*1 

=  n" 

+  cos  A" 

y2 

=  1 

+  n"  cos  a" 

6" 

=Y7 

-  e"2 

r 

1_ 

a 17 

-  s 

a 

=  (1 

-  n’’V/2"c 

CD 

=  aerodynamic  drag 

S4 


A  =  effective  cross-sectional  area  of  satellite 

m  =  mass  of  satellite 

B  =  i  Cn  -  .  ,qo  -  s, t 
°  2  D  m  PQ  q,< 


B  =  B_  (a”  -  s1T 
o 

0"  =  COS  I” 

n"  =  V  y  (a")  =  "mean"  mean  motion 

1  2 

k2  =  2  J2  R  where  R  is  the  Earth's  average  equatorial 
radius  and  3^  is  the  second  zonal  coefficient  in 
the  geopotential 

A3  0  "  '  J3  where  is  the  third  zonal  coefficient 

in  the  geopotential 
7i  4 

k4  =  ’  8  J4  R  where  J4  is  the  forth  zonal  coefficient 
in  the  geopotential 

C 

A$  0  =  -  Jj  R  where  Jj-  is  the  fifth  zonal  coefficient 
in  the  geopotential 


I . 


3  ,k 


1 

W 


I 


2  7T 


cos  jx 


(1  - n"  cos  x) 


dx 


-  Cn")J  £ 

(k  -  1)1  1 


where  ^  =  Y  1  -  n 


=  </> 


-1 
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f2  =  3^.  2  +  ^  3 

f3  =  (j2  -l)'!'  3  +  3jip  4  +  3ijj  5 

f4  =  (j3  -4JH'4  +  (6j2  -9)*'5  +  3  Sjip" 6  +  ISf7 

f5  =  (j4  -10j2  +  9)>jr5  +  (10j3  -  5 S j ) ip " ^ 

+  (45 j 2  -  90)ip"7  +  lOSjip*8  +  105ip"9 
f6  =  (j5  -  20  j  3  +  64j)ip"6  +  (15j  4  -  195  j  2  +  225)ip 
+  (105 j  3  -  735j)tp"8  +  (420 j 2  -  lOSO)^'9 
+  945ji//10  +  94  Sip  1  ^ 

f?  =  (j6  -  35  j 4  +  2  59  j  2  -  225)ip  * 7  +  (21j5-  52  5  j  3 
+  2079 j  ) 8  +  (2 10 j  4  -  3360 j  2  +  4725)ip'9 

+  (1260  j  3  -  10 , 7 1 0  j  )ip"10  +  ( 4  7  2  5  j  2  -  14,17S)i|> 

+  10, 395jip  12  +  10,395ip'13 
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